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Abstract 

Let p(x 1 ,...,x n ) = J2( ri ,...,r n )ei„,„ a (r 1 ,...,r n ) Ui<i< n x 4 be homogeneous polynomial of 
degree n in n real variables with integer nonnegative coefficients. The support of such 
polynomial p(x\, x n ) is defined as supp(p) = {(ri,...,r„) S I n ,n '■ a (ri,...,r n ) 7^ 0} . The 
convex hull CO(supp(p)) of supp(p) is called the Newton polytope of p . We study the 
following decision problems , which are far-reaching generalizations of the classical perfect 
matching problem : 

• Problem 1 . Consider a homogeneous polynomial p{x\, ...,!„) of degree n in n real 
variables with nonnegative integer coefficients given as a black box (oracle ) . Is it 
true that (1, 1, .., 1) G supp(p) ? 

• Problem 2 . Consider a homogeneous polynomial p(xi, ■ ■■,x n ) of degree n in n real 
variables with nonnegative integer coefficients given as a black box (oracle ) . Is it 
true that (1, 1, .., 1) G CO(supp(p)) ? 

We prove that for hyperbolic polynomials these two problems are equivalent and can be 
solved by deterministic polynomial-time oracle algorithms . This result is based on a "hy- 
perbolic" generalization of the Rado theorem . We also present combinatorial and algebraic 
applications of this " hyperbolic" generalization of the Rado theorem (prove that the support 
supp{p) of P-hyperbolic polynomial p is an intersection of some Integral Polymatroid with 
the hyperplane {(r±, r n ) : J2i<i<n r i = n } ) an d pose some open problems. 



1 Introduction and motivating examples 

The layout of the paper : 

We introduce the main topics and motivations in Section 1 . In Section 1.1 we present a naive 
algorithm to solve Problem 1 in the general case . We show in Appendix D that this algorithm 
is , in a sense , optimal . 

(Incidentally (or not) , the situation here is very similar with the optimality of the square root 
in the famous quantum Grover's search algorithm .) 

In Section 1.2 we remind the basic properties of hyperbolic polynomials used in this paper . 

In Section 2 we state a hyperbolic analogue of the Rado theorem (Theorem 2.2) , which is the 
main mathematical result of the paper . Theorem 2.2 sheds more light on the algebraic- geometric 
nature of such fundamental combinatorial results as Hall 's and Rado 's theorems . 
In Section 2.1 we define and study doubly-stochastic polynomials (an useful generalization of 
standard doubly-stochastic matrices ). We also state there a hyperbolic analogue of the van der 
Waerden conjecture . 
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In Section 3 we introduce and analyse the ellipsoid algorithm which solves Problem 1 and Prob- 
lem 2 on the class of S -hyperbolic polynomials . The essence of the results in Section 3 is that 
once Hall's like conditions (the exponential number of them) are proved to be necessary and 
sufficient , they can be checked by a polynomial time deterministic oracle algorithms . The 
algorithm , which we use , is based not on the linear programming but on some rather nonlinear 
convex programs similar to considered in \23)[ , 1211/ , \22j . 

In section 4 we introduce and analyse another algorithm , which is a "polynomial" general- 
ization of the Sinkhorn Scaling . 

In Section 5 we use Theorem 2.2 to get a more refine (polymatroidal) properties of the supports 
of P -hyperbolic polynomials and explain how our results generalize the main result from 

In Section 6 we pose some open problems and share our enthusiasm about the topic . 
The proofs of the main results are presented in Appendices A,B,C,D . 

Let ...,x n ) = E(ri,...,r„)eJ n ,n a (n,-,0 Ui<i<n x ? De homogeneous polynomial of degree 
n in n real variables. Here Ik } n stands for the set of vectors r = (Vi, r^) with nonnegative inte- 
ger components and J2i<i<k r i = n - ^ n this paper we primarily study homogeneous polynomials 
with nonnegative integer coefficients . 

Definition 1.1: The support of the polynomial p(xi, ...,x n ) as above is defined as supp(p) = 
{(ri, r n ) G I n ,n '■ 0'(n,...,r n ) 7^ 0} • The convex hull CO(supp(p)) of supp(p) is called the 
Newton polytope of p . I 

We will study the following decision problems : 

• Problem 1 . Consider a homogeneous polynomial p(x±, ...,x n ) of degree n in n real 
variables with nonnegative integer coefficients given as a black box (oracle ) . Is it true 
that (1, 1, .., 1) 6 supp(p) ? 

An equivalent question is : Is it true that q x ® q Xi -p(xi, ...,£„) ^ ? 

• Problem 2 . Consider a homogeneous polynomial p(xi, x n ) of degree n in n real 
variables with nonnegative integer coefficients given as a black box (oracle ) . Is it true 
that (1, 1, .., 1) 6 CO(supp{p)) ? 

Our goal is solve these decision problems using deterministic polynomial-time oracle algorithms 
, i.e. algorithms which evaluate the given polynomial p{.) at a number of rational vectors 
(qi, q n ) which is polynomial in n and log(p(l, 1, .., 1)); these rational vectors (q%, q n ) are 
supposed to have bit-wise complexity which is polynomial in n and log(p(l, 1, .., 1)) ; and the 
additional auxilary arithmetic computations also take a polynomial number of steps in n and 

logfrCM,..,!)) • 

The next example explains some (well known ) origins of the both problems . 
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Example 1.2: Consider first the following homogeneous polynomial from [2S] : p(%i, ■■■>x n ) = 
tr{{D{x)A) n ) , where D{x) is a n x n diagonal matrix Diag(x\, x n ) ; and A is n x n ma- 
trix with (0, 1) entries , i.e. A is an adjacency matrix of some directed graph T . Clearly 
, this polynomial p(x\, x n ) has nonnegative integer coefficients . It was proved in [2Hj 
that ^ q x ^"q x tr((D(x)A) n ) is equal to the number of Hamiltonian circuits in the graph T 
. Notice that the polynomial tr(D(x)A) n can be evaluated in 0(n 3 log(n)) arithmetic oper- 
ations and (1,1,.. .,1) 6 supp(p) iff there exists a Hamiltonian circuit in the graph T. Also 
log(p(l, 1, .., 1)) < nlog(n). Therefore , unless P = NP , there is no hope to design determin- 
istic polynomial-time oracle algorithm solving Problem 1 in this case . (The author is indebted 
to A.Barvinok for pointing out this polynomial . ) 

Consider , with the same adjacency matrix A , another homogeneous polynomial 
Mul(x!,...,X n ) =Ul<i< n T,l<i<n A ( i ,j) x j • Then 
1J £^-Mul(x 1 ,...,x n ) = Per(A) . 

Therefore for the multilinear polynomial Mul(x±, ...,x n ) Problem 1 is a "black box" analogue 
of checking the existence of the perfect bipartite matching . 



Next consider the following class of determinantal polynomials : 

q(xi, ...,x n ) = det( A { Xi), 

l<i<n 

where A = (A\, ...,A n ) is a n-tuple of positive semidefinite n x n hermitian matrices , i.e. 
Ai y , with integer entries . Recall that the mixed discriminant 

Z) ( A ) = ^ ^ det ( E A ^)- 

da x ...da n ^ n 

(If the matrices Ai above are diagonal , i.e. Ai = Diagipi^, bi,n), 1 < i < n , 

then the mixed discriminant is reduced to the permanent : D(A) = Per(B), B = {pij;l < 

i,j <n ). 



It is well known (see , for instance , [22] ) that a determinantal polynomial q(.) can be 
represented as 

q(x u ..,x n )= ]T J] x ? D ( A r) u . 1 (!) 

r£l ni „l<i<n lll<i<n r «- 

where a n-tuple A r of square matrices consists of r, L copies of Ai, 1 < i < k . One of the 
equivalent formulations [25 °f the classical Rado theorem states that D (An ^ > iff 

Rank(J2 A i) > \ S \ for a11 S c i 1 ^ 2 > •••»"} ( 2 ) 

ieS 

(The diagonal case is the famous Hall's theorem on the perfect bipartite matchings .) 
The Rado theorem is just a particular case of famous Edmonds theorem on the rank of inter- 
section of two matroids . Therefore , given a n-tuple A = (A\, A n ) of positive semidefinite 
n x n hermitian matrices , one can decide in deterministic polynomial time if D(A) > . We 
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will explain below that this decision problem can be solved by a deterministic polynomial-time 
oracle algorithm . I.e. we only use some values of det(Xa<i< n M^i) without reconstructing the 
actual tuple A = (Ai, A n ) . 

The natural question , in our opinion , is which algebraic-geometric properties 
make the class of determinantal polynomials "easy" and the class of Barvinok's 
polynomials tr(D(x)A) n "hard" . This paper suggests one answer to the question. 

One important corollary of the Rado conditions (2) is that 

supp(q) = CO{supp(q)) n /„,„. (3) 
I.e. if integer vectors r, r(l), r(2), r(k) G I(n,n) and 

r = a(i)r(i),a(i) > 0, 1 < i < k; a(z), 

l<i<k \<i<k 

and D(A r (j)) > 0, 1 < i < k then also D(A r ) > . Notice that in this case Problem 1 and 

Problem 2 are equivalent . 

We can rewrite Rado conditions (2) as follows : 

max Vr; > |5| for all Sc {l,2,...,n} (4) 

r£supp(q) i( _ s 

Putting things together we get the following Fact . 

Fact 1.3: The following properties of determinantal polynomial q(x±, x n ) = det(J2i<i< n -^i x i) 
with n x n hermitian matrices Ai >z 0, 1 < i < n are equivalent . 

1. (1,1, ..,1) i supp{q). 

2. (1,1,.,1) <£CO(supp(q)). 

3. There exists nonempty S C {1,2, ...,n} such that 

r i s i< Sj = IS 1 ! for all(ri, ...,r n ) G supp(^), (5) 

l<i<n l<i<n 

where (s\, ...,s n ) is a characteristic function of the subset S , i.e. Sj = 1 if i G S , and 
Si = otherwise . 

Notice that if (5) holds then the distance dist(e, CO{supp{q))) from the vector e = (1, 1) 

to the Newton polytope 

CO{supp{q)) is at least ^ \ s \ { n-\s\) ^ ^ ■ 

I 
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We will show that for any class of polynomials satisfying Fact 1.3 there exists a deterministic 
polynomial-time oracle algorithm solving both Problem 1 and Problem 2 , which are , of course 
, equivalent in this case . Our algorithm is based on the reduction to some convex programming 
problem and the consequent use of the Ellipsoids method . 

The next fact about determinantal polynomials , namely their hyperbolicity , is "responsible" 
for Fact 1.3 . 

Fact 1.4: Consider a determinantal polynomial q((x±, x n ) = det(J2i<i<n A % x i) with Ai >z 
0, 1 < i < n . Assume that q is not identically zero , i.e. that B =: J2i<i< n A i >- (the sum is 
strictly positive definite ). For a real vector {x\, ...,x n ) G R n consider the following polynomial 
equation of degree n in one variable : 

P(t) = q(xi -t,x 2 ~t,...,x n - t) = det( A i x i-t A i)=°- ( 6 ) 

l<i<n l<i<n 

Equation (6) has n real roots roots counting the multiplicities ; if the real vector (xi,...,x n ) G R n 
has nonnegative entries then all roots of (6) are nonnegative real numbers . | 

The main result of this paper that this hyperbolicity , which we will describe formally in 
Section 1.1 , is sufficient for Fact 1.3 ; i.e. Fact 1.4 implies Fact 1.3 . I 

1.1 "Naive" algorithms 

One possible " naive" algorithm to solve Problem 1 is just to compute q x ^ q x Pi^ii •••^x n ) 
. Recall that the number of coefficients of a homogeneous polynomial of degree n in n real 
variables is equal to « 2 2n . We can compute all the coefficients of p(x±, x n ) via 

standard multidimensional interpolation , but this interpolation will need ^^-\)\ ~ ^ 2ra oracle 
calls . There is an algorithm which computes dx ® n dx p{x\, —,x n ) using only 2" _1 oracle calls : 

^£^-p( Xl ,...,x n )= (7) 

2 _n+1 E 6ie {-l, + l},2<K» l J ) ( 1 , bn) I\2<i<n b i- 

This formula is , in a sense , optimal . I.e. there exists a nearly matching lower bound . The 
corresponding result and connections to computations of the permanent are presented in Ap- 
pendix D . 

We will explain below that if p is a homogeneous polynomial of degree n in n real variables 
with nonnegative integer coefficients then (1, 1, .., 1) G CO(supp(p)) iff p(x±, x n ) > rii<i<n x « 
for all vectors ...,x ra ) with positive real coordinates . Therefore Problem 2 is equivalent to 
checking if the polynomial P(y u ...,y n ) = p(l + yf, 1 + yl) - Ui<i< n 1 + Vi is nonnegative 
on R n . 
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1.2 Hyperbolic polynomials 



The following concept of hyperbolic polynomials was originated in the theory of partial differ- 
ential equations [TH], [5] ,^U| . 

A homogeneous polynomial p(x),x G R m of degree n in m real varibles is called hyperbolic in 
the direction e G R m (or e- hyperbolic) if for any x G i? m the polynomial p(x — Ae) in the one 
variable A has exactly n real roots counting their multiplicities. We assume in this paper that 
p(e) >0. Denote an ordered vector of roots of p(x—Xe) as X(x) = (X\(x) > X2(x) > ...X n (x)). It 
is well known that the product of roots is equal to p(x). Call x G R m e-positive (e-nonnegative) 
if X n (x) > (A n (x) > 0). The fundamental result JH] hi the theory of hyperbolic polynomi- 
als states that the set of e-nonnegative vectors is a closed convex cone. A /c-tuple of vectors 
(xi, ...Xk) is called e-positive (e-nonnegative) if Xi, 1 < i < k are e-positive (e-nonnegative). We 
denote the closed convex cone of e-nonnegative vectors as N e (p), and the open convex cone of 
e-positive vectors as C e {p). 

Recent interest in the hyperbolic polynomials got sparked by the discovery that 
\og(p(x)) is a self- concordant barrier for the opened convex cone C e (p) and therefore the pow- 
erful mashinery of interior-point methods can be applied . It is an important open problem 
whether this cone C e (p) has a semi-definite representation . 

It has been shown in ^H] (see also |26|) that an e- hyperbolic polynomial p is also d- 
hyperbolic for all e-positive vectors d G C e {p). 

Let us fix n real vectors x% G R m , 1 < i < n and define the following homogeneous polynomial: 

Px u ..,x n {ai,...,a n ) =p( a i x i) ( 8 ) 

\<i<n 

Following [201 j we define the p-mixed form of an n- vector tuple X = (x%, -,x n ) as 

d n 

M P (X) =: M p (xi, ..,x n ) = — p( 22 diXi) (9) 

1 "' n l<i<n 

Equivalently, the p-mixed form M p (x\, ..,x n ) can be defined by the polarization (see [26 ) : 

Mp( Xl ,..,x n ) =2- n P(J2 b * x i) II b i ( 10 ) 

b«e{-l,+l},l<i<ra l<«<n l<i<n 

Associate with any vector r = (n,...,r n ) G l n ^ n & n n-tuple of m-dimensional vectors X r 
consisting of rj copies of X{(1 < i < n). It follows from the Taylor's formula that 

P Xl ,.., Xn ( ai ,...,a n )= II «?M p (X r ) 1 (11) 

re/„,„l<i<n lll<i<n r i- 

For an e-nonnegative tuple X = (xi, define its capacity as: 

Cap(X) = inf Pxi,..,x n ( a U «n) 

a ' >0 'T\i<i<n a > =1 
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Probably the best known example of a hyperbolic polynomial comes from the hyperbolic 
geometry : 

P(a ,...,a k ) = al - ^ a- (13) 

i<j<fc 

This polynomial is hyperbolic in the direction (1,0,0, ...,0). Another "popular" hyperbolic 
polynomial is det(X) restricted on a linear real space of hermitian nx n matrices . In this case 
mixed forms are just mixed discriminants , hyperbolic direction is the identity matrix / , the 
corresponding closed convex cone of J-nonnegative vectors coincides with a closed convex cone 
of positive semidefinite matrices . 

Less known , but very interesting , hyperbolic polynomial is the Moore determinant Det(M\(Y) 
restricted on a linear real space of hermitian quaternionic n x n matrices . The Moore deter- 
minant is , essentially , the Pfafhan (see the corresponding definitions and the theory in a very 
readable paper [HE] ) • 

We use in this paper the following class of hyperbolic in the direction (1, 1, 1) polynomials 
of degree k : 

Q(ax,...,a k ) = M p (J2i<i<k a i x i, -,T,i<i<k a i x ii x k+l, -,%n), where p is a e-hyperbolic poly- 
nomial of degree n > k , (x\, ..,x n ) is e-nonnegative tuple , and the p-mixed form 

Mp(J2l<i<k x ii J2l<i<k x ii x k+l, •••> X n ) > 0. 

We make a substantial use of the following very recent result [22j > w hich is a rather direct 
corollary of , [H7j . 

Theorem 1.5: Consider a homogeneous polynomial p(y 1,1/2, 2/3)) of degree n in 3 real variables 
which is hyperbolic in the direction (0,0,1). Assume that p(0, 0, 1) = 1 . Then there exists two 
n x n real symmetric matrices A, B such that 

p(yi,V2, 2/3)) = det(yiA + y 2 B + y 3 I). 

It has been shown in ^H] that most of known facts , and some opened problems as well , about 
hyperbolic polynomials follow from Theorem 1.5 . 

2 A hyperbolic analogue of the Rado theorem 

Definition 2.1: Consider a homogeneous polynomial p{x),x G R m of degree n in m real 
variables which is hyperbolic in the direction e. Denote an ordered vector of roots of p(x — Xe) 
as X(x) = (Ai(x) > A2(x) > ...X n (x)) . We define the p-rank of x £ R m in direction e as 
Rank p {x) = \{i : Xi(x) 7^ 0}|. It follows from Theorem 1.5 that the p-rank of x G R m in any 
direction d G C e is equal to the p-rank of x G R m in direction e , which we call the p-rank of 
x G R m . I 

Consider the following polynomial in one variable D(t) = p(td + x) = J2o<i< n c $ ''■ ^ follows 
from the identity (11) that 
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c n = M p (d,..,d)(nl)- 1 =p(d), (14) 
c n _! = M p (x,d, ., d)(l\(n - l)!)- 1 , .., 
Co = M p (x, .., x){n\)~ l = p{x). 

Let (A^(x) > \'^'\x) > ... > A^^x)) be the (real) roots of x in the e-positive direction d, i.e. 
the roots of the equation p(td — x) = . Define (canonical symmetric functions) : 

Sk,d( x ) = Xi 1 (x)X i2 (x)...X ik (x). 

I<il<i2<...<ik<n 

Then Sk,d{ x ) = ■ Clearly if x is e-nonnegative then for any e-positive vector d the p-rank 
Rankp(x) = max{&; : Sfc i( i(x) > 0} . The next theorem , which we prove in Appendix A , is 
the main mathematical result of this paper . Our main tool is Theorem 1.5 , which facilites a 
rather easy induction . We also use a particularly easy case of the Rado theorem (see Remark 
A. 5 for the details ) . 

Theorem 2.2: Consider a homogeneous polynomial p(x),x £ R m of degree n inm real variables 
which is hyperbolic in the direction e , p(e) > 0. Let X = (x±, ...x n ),Xi € R m be e-nonnegative 
n-tuple of m- dimensional vectors , i.e. Xj,l < i < n are e-nonnegative . 

Then the p-mixed form M p (X) =: M p (x±, ..,x n ) is positive iff the following generalized Rado 
conditions hold : 

RankpC^Xi) > \S\ for all S C {1, 2, n}. (15) 

ies 

Definition 2.3: Call a homogeneous polynomial p(a),a G R n of degree n in n real variables 
P- hyperbolic if it is hyperbolic in direction e = (l,l,...l) (vector of all ones) , p(e) > and all 
the canonical orts e^, 1 < i < n (rows of the identity matrix / ) are e-nonnegative . In other 
words ,a homogeneous polynomial p(ct), a £ R n of degree n in n real variables is P-hyperbolic if 
it is e-hyperbolic and its closed cone of e-nonnegative vectors contains the nonnegative orthant 
RJ] = {(xi, x n ) : Xi > 0, 1 < i < n} . It follows from |26| that the coefficients ofP-hyperbolic 
polynomials are nonnegative real numbers . 

(Notice that the class of P-hyperbolic polynomials coincides with the class of polynomials 
) = P(J2i<i<n a i x i) > where p is e -hyperbolic polynomial of degree n in m 
real variables , a n-tuple (xi, ..,x n ) of m-dimensional real vectors is e-nonnegative and J2i<i< n x i 
is e-positive . ) 

Call a homogeneous polynomial q(a),a G R n of degree n in n real variables with nonnegative 
coefficients S-hyperbolic if there exists a P-hyperbolic polynomial p such that supp(p) = supp(q) 
I 

(One natural class of S'-hyperbolic polynomials ,not all of them P-hyperbolic , is Vol(a\Xi + 
...a n X n ) , where Xi are convex compact subsets of R n . See the corresponding not P-hyperbolic 
example in [^S] .) 
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Corollary 2.4: Let q(a),a G R n be S '-hyperbolic polynomial of degree n . 
Then CO{supp(q)) l~l I n>n = supp(q) . 

Proof: It is enough to prove the corollary for P-hyperbolic polynomials. I.e. suppose that 
q(a±, a n ) = p(J2i<i<n a % x i) ■> where p is e -hyperbolic polynomial of degree n in m real 
variables , a n-tuple (xi,..,x n ) of m-dimensional real vectors is e-nonnegative and J2i<i<n x i 
is e-positive . Then r = (r±,r2, ■■■,r n ) G supp{q) iff the p-mixed form M p (X r ) > , where the 
n-tuple X r consists of ri copies of Xi, 1 < i < n. Let r^ = (rf \ ...,ri°' ) ) G CO{supp{q)). I.e. 
there exist r^) G supp(q), 1 < j <n such that r' ' = X)i<j<n a j r ^ an< ^ a i — 0, X)i<j<n a i = 1 

Let rtf) = (ri i} , ri i} ), < j < n . As G supp(q), l<j<n thus M p (X r0) ) > 0, 1 < j < n 
. It follows from Theorem 2.2 (only if part ) that 

Rank p {^2 %i) > E r ? } for a11 5 c i 1 ' 2 ' -' n ^ 1 - 3 - n - 

i&S i&S 

Therefore 

i?an£; p £»>E E S C {1,2, ..,n}. 

i£S i£Sl<j<n i£S 

Using the "if part of Theorem 2.2 we get that M p (X r ( () )) > and thus r(°) G supp(q) . I 

Corollary 2.5: Let q(x),x G 6e S-hyperbolic polynomial of degree n . Then the following 
conditions are equivalent 

1. e G CO(supp(q)) . 

2. e G supp{q) , i.e. da ^ dan q{x) > . 

3. Cap{p) =: inf tt ai=1 <z(ai, .... a n ) > 0. 

^. For all e > there exists a vector (ai, ...,a n ) i/rti/i positive entries such that the following 
inequality holds : 

da ; ^-l 2 <e. (16) 

!<*<n 9("l.-.«n) 

5. There exists a vector (a±, a n ) with positive entries such that the following inequality 
holds : 

a < ^ g (ai,...,a w ) ^ 1 
i^<„ ^(«i,-,a„) -n" 

6. For a// subsets S C {1,2, ...,n} the following inequality holds : 

E r * - \ S \ f° r al1 ( r i'---> r n) G supp{q). (18) 
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(We sketch a proof in Appendix C . ) 



The following result , which we prove in Appendix B , is a "polynomial" generalization of 
Lemma 4.2 in (20] . 

Proposition 2.6: The condition (17) implies the condition (18) for all homogeneous polyno- 
mial q{x),x £ R n of degree n in n real variables with nonnegative coefficients . 

2.1 Doubly-stochastic polynomials 

Inequalities (16), (17) above suggest the following definition . 

Definition 2.7: A homogeneous polynomial q{x\, x n ) of degree n in n variables is called 
doubly-stochastic if its coefficients are nonnegative real numbers and ^-g(l, 1, 1) = 1 for 
all 1 < % < n. The doubly-stochastic defect of the polynomial q is defined as DS(q) = 
El<i<n(^<?(M, •••,!)- I) 2 ■ 



Lemma 2.8: 



1. A homogeneous polynomial q(x±, x n ) of degree n in n variables with nonnegative real 
coefficients is doubly- stochastic iff q(l, 1, 1) = 1 and q(x±, ...,x n ) > Yli<i< n Xi for all 
real vectors (xi, x n ) £ R n with positive coordinates (in other words if q(l, 1, 1) = 1 
and Cap(q) = 1 ). 

2. A homogeneous polynomial q(x\, x n ) of degree n in n variables is P -hyperbolic and 
doubly-stochastic iff q(l, 1, 1) = 1 and \q(z±, z n )\ > JJi <i<n Re{z,{) for all complex 
vectors vectors (zi,...,z n ) G C n with positive real parts . 

3. If a sequence qi of homogeneous polynomials of degree n in n variables with nonnegative 
real coefficients converges to a doubly- stochastic polynomial then lim^oo Cap(qi) = 1. 

4- The capacity Cap(q) is a continuous functional (but not even Lipshitz) on a convex closed 
cone of homogeneous polynomials of degree n in n variables with nonnegative real coeffi- 
cients . 

5. If q is homogeneous polynomial of degree n inn variables with nonnegative real coefficients 
then Cap(q) > iff the exists a sequence Xj = (xij, ...,x n j) of vectors with positive real 
coordinates such that 

■■■■,x n ,j) ^ 2 ^ 

i<i<n —,x n ,j) 

And in this case Cap(q) = lim^oo ^ q<yX ^ — . 

lh<i< n Xi >3 
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Example 2.9: A multilinear polynomial 

q(xi, ...,x n ) = Y\i<i< n J2i<j<n a (h j) x j is doubly-stochastic and P-hyperbolic iff 

the square matrix B = — a ( t <^) ■ \ < i,j < n} is doubly stochastic in the standard 

meaning . | 

The next theorem is another Corollary of Theorem 2.2 ; its main point is in introducing a 
hyperbolic analog of Van der Waerden conjecture . 

Theorem 2.10: 

1. Consider the set PHDS(n) of all P-hyperbolic doubly-stochastic homogeneous polynomials 
q(xi, x n ) of degree n in n variables. The set PHDS(n) is a compact and the following 
inequality holds 

inf ge PHDS(n) dxi...dx n q ( Xl >-> Xn ) = 
™n- qePHDS{n) dxi d " dXn q(xi,...,x n ) =: VdW(n) > 

2. For any P-hyperbolic homogeneous polynomials q(x\, ...,x n ) of degree n in n variables the 
following inequalities hold 

VdW{n) < 9 Xl ... 9x ^--^n) L 
Cap{q) 

Conjecture 2.11: Hyperbolic Van der Waerden conjecture 

VdW(n) = — ? 



The next result , a direct corollary of Lemma 2.10 in ^S] , is a generalization of Proposition 
4.2 in 

Lemma 2.12: Let q be P-hyperbolic homogeneous polynomial of degree n inn variables . Then 
the following inequalities hold for all vectors (x±, ...,x n ) with positive real coordinates. 



1. 

d 



lll<i<n x i Ki< n yV x l) x n) 



2. //log( ^ xl '-' x " ) ) - \og(Cap(q)) < ^ with < e < 1 then 

1 ll<i<n X% 



\^±r. ^ e (2°) 



n 



Example 2.13: Consider the following homogeneous polynomial of degree n in n variables : 
p(x 1 , ...,x n ) = J2i<i<n x i- Tli en Cap(p) = n and 

P(»l,-,^) n £- P (x u ...,x n ) _ 

= ""( j! 1 ^': )"- 1 < n = Cap{p). 

The moral of this example is that the inequality (19) does not hold for all homogeneous poly- 
nomials with nonnegative coefficients , this inequality is a nontrivial necessary condition for the 
P-hyperbolicity . It is interesting to notice that the inequality (19) implies the determinantal 
Hadamard inequality . I 

Remark 2.14: Let 

p(xi,...,x n )= a {n,..,r n ) II x ? 

(ri,...,r n )e/n,n l<i<71 

be homogeneous polynomial of degree n in n real variables. Assume that its coefficients are 
nonnegative and sum to one , i.e. that p(l, 1, 1) = 1. Associate with this polynomial a 
random integer vector 

Z p = (zi,...,z n ) G l n ,n ■ Prob{Z p = (n,...,r n )} = a( rir .. )f . n ). 

Then 

d d 

Therefore p is doubly-stochastic iff E(Z p ) = (1, 1, 1) ; there exists a doubly-stochastic poly- 
nomial q such that supp(q) G supp(p) iff (1, 1, 1) G CO(supp(p)). 

It follows from Corollary 2.5 that if p is doubly-stochastic and S- hyperbolic then Prob{Z p = 
E(Z p )} > . And the hyperbolic van der Waerden conjecture can be reformulated as : 
If p is doubly-stochastic and P -hyperbolic then 
Prob{\\Z p = E(Z P )\\ <V2}>£ . 

Perhaps some kind of the measure concetration is present here ? | 



Remark 2.15: The problem to find out a positive real solution of the inequality (20) is a 
far reaching generalization of scaling of matrices with nonnegative entries (the corresponding 
polynomials are multilinear )[10] , ESI an d scaling of tuples of PSD matrices (the corresponding 
polynomials are determinantal ) [2J| > |22j • Part 2 of Lemma 2.12 allows to generalize results 
of |23j , |21| , |22| to P-hyperbolic polynomials , even in the black-box setting . Can it be done 
for all homogeneous polynomials with , say , integer nonnegative coefficients ? I 



3 The ellipsoid algorithm 



Consider a homogeneous polynomial q(x),x G R n of degree n in n real variables with nonneg- 
ative integer coefficients . Associate with such q the following convex functional 

f(y h ...,y n )=log(q(ey\eV*,...,e^)). 
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Proposition 3.1: The following conditions are equivalent 



1. e= (1,1,..,1) G CO(supp(q)) . 

2. mf yi+ ... +yn= of(yi,...,y n ) > 0. 

3. Ife= (1, 1, .., 1) i CO(supp(q)) then 
inf ai+ ... +2/n=0 f(yi, ...,y n ) = -oo. 

Let dist(e, CO{supp{q))) = A" 1 > and Q = log(g(e)) . De/ine 7 = (Q + 1)A . T/ien 

illf j_ _l rwi I2_l ,1 12^ /(2/l>->2/«) = 
yi+...+j/„=0,(|yi| i +... + |j/„| i )^<7 

min 2/1+ ... +2/n=0! | 2/1 |2 + .„ + | 2/n |2< 7 /(yi J ...,?/ n ) < -1. 



Proof: Our proof is a straigthforward application of the concavity of the logarithm on the 
positive semi-axis and of the Hanh-Banach separation theorem . It will be included in the full 
version . | 

Proposition 3.1 suggests the following natural approach to solve Problem 2 , i.e. to decide 
whether e = (1, 1, .., 1) G CO{supp(q)) or not : 

find nrim^ ...+y n =o,\ yi p+ ...+\y n \ 2 <-y fiv^i ■■■i Vn) with absolute accuracy | . If the resulting value 
is greater than or equal —3 then e = (1, 1, .., 1) G CO(supp(q)) ; if the resulting value is less 
than or equal — | then e = (1, 1, .., 1) ^ CO(supp(q)) . And , of course , it is natural to use the 
ellipsoid method . Our main tool is the following property of the ellipsoid algorithm [32]: For a 
prescribed accuracy S > 0, it finds a 5-minimizer of a differentiable convex function / in a ball 
B, that is a point xg £ B with f(x$) < min^ / + 5, in no more than 

O In In I II, (Var B (/) = max/ - nun/) 

iterations. Each iteration requires a single computation of the value and of the gradient of / at 
a given point, plus 0(n 2 ) elementary operations to run the algorithm itself. In our case, this is 
easily seen to cost at most 0(n 2 ) oracle calls and 0(n) elementary arithmetic operations . 
We have the n — 1-dimensional ball B 1 = {(yi, ■■■,y n ) ■ Vi + ••• + Vn = 0, \y\\ 2 + ... + \y n \ 2 < 7}- 
A straigthforward computations show that 

Var B (f) < log(g(l, 1, .., l)e^) - log(g(l, 1, .., l)e"^) < 2 7 n, 

giving that 0(n 2 (ln(n) +ln(7)) iterations of the ellipsoid method needed to solve Problem 2 , it 
amounts to 0(n 4 (ln(n) +ln(7)) oracle calls . The quantity 0(n 4 (ln(n) +ln( 7 )) is polynomial in 
n even if 7 is exponentially large (dist(e, CO(supp(q))) is exponentially small ). The problem 
is that if 7 is exponentially large ( which can happen ) then we need to call oracles on inputs 
with exponential bit-size . 

Putting things together , we get the following conclusion : 

If it is promised that either e = (1, 1, .., 1) G CO(supp{q)) or dist{e, CO(supp(q))) > poly{n)~ l 
for some fixed polynomial poly{n) then Problem 2 can be solved by a deterministic polynomial- 
time oracle algorithm based on the ellipsoid method . 

And at this point we can say nothing about Problem 1 , i.e. deciding whether e = (1, 1, .., 1) G 
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supp(q) or not . Corollary 2.5 says that if q is .S-hyperbolic polynomial then Problem 1 and 
Problem 2 are equivalent ; moreover if e = (1,1, ..,1) ^ supp(q) then here exists nonempty 
S C {1, 2, n} such that 

Y r i s i < Y Si = \S\ iov a\l(ri,...,r n ) e supp(q), (21) 

l<i<n l<i<n 

, where (si, s n ) is a characteristic function of the subset S , i.e. Sj = 1 if i G S , and s« = 
otherwise . 

Notice that if (21) holds then the distance dist(e, CO{supp{q ))) from the vector e = (1,...,1) 
to the Newton polytope CO(supp(q)) is at least \ s \ (n- \ g \j — "^k ' Thus we have the next 
theorem . 



Theorem 3.2: Problem 1 and Problem 2 are equivalent for S-hyperbolic polynomials . 

There exists a deterministic polynomial-time oracle algorithm solving Problem 1 for a given 

S-hyperbolic polynomial q(ai, ■■■,ct n ) with integer coefficients . 

It requires 0(n 4 (ln(n) + ln(ln(g(l, 1, 1))) oracle calls and it bit-wise complexity (which is 
roughly the radius of the ball ) is 0(n? ln(g(l, 1, 1))) . 

4 Hyperbolic Sinkhorn scaling 

We will discuss briefly in this section another method , which is essentially a large step version 
of gradient descent . 

Definition 4.1: Consider an e-nonnegative tuple X = (xi,..,x n ) such that the sum of its 
components S'(X) = d = J2i<i<k x i 1S e-positive. Define tr^(x) as a sum of roots of the 
univariate polynomial equation p(x — td) = 0. 

Define the following map (Hyperbolic Sinkhorn Scaling) acting on such tuples: 

HS(X) = Y = ( — ^ — -, — ^r) 

V ^ y tr d { Xi y 'tr d {x n )> 

Hyperbolic Sinkhorn Iteration (HSI) is the following recursive procedure: 

X J+ i = HS(X.j),j > 0, Xo is an e-nonnegative tuple with 
J2l<i<k Xi £ C e . 

We also define the doubly-stochastic defect of e-nonnegative tuples with e-positive sums as 

DS(X) = ]T (tr d ( Xi ) - l) 2 ; ]T x t = d G C e 

l<i<k l<i<k 

I 

We can define the map HS(.) directly in terms of the P-hyperbolic polynomial 
Q(ai,...,a n ) = P Xl ,..,x n (ai,-,a n ) = p( aiXi )- 

l<i<n 
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Indeed, if J2 



l<i<n u -i- L l 



ctiXi = d E C e then 



tr d {a.iXi) 



aiQ^-Q{ai,...,a n ) 



(22) 



Q(ai,...,a n ) 



This gives the following way to redefine the map HS(X.) : 



HS(ai, ...,a n ) = ( 



...,a n ) 



Q(a 1 ,...,a n ) 



^-Q(ai,...,a n ) 



^Q{a u ...,a n ) 



for aj > 0, 1 < % < n. 

And correspondingly the doubly-stochastic defect of (ai, ■■■,ot n ) is equal to 




otiQ^-Q(ai,...,a n ) 
Q(ai, ...,a n ) 



the same as the left side of (17 ) . Notice that J2i<i<n ^ r d{xi) = n by the Euler's identity . 

Example 4.2: Consider the following hyperbolic polynomial in n variables: p{z\, z n ) = 
Yli<i<n z i- ^ is e- hyperbolic for e = (1,1,.. .,1). And N e is a nonnegative orthant, C e is a 
positive orthant. An e-nonnegative tuple X = (xi, .., x n ) can be represented by an n x n matrix 
Ax. with nonnegative entries: the ith column of A is a vector Xi G i? n . If Z = (z%, z n ) € .R n 
and d = (d\, ...,d n ) £ R n ; Z{ > 0, 1 < i < n, then trd(Z) = J2i<i<n Recall that for a square 
matrix A = {ajj : 1 < i,j < N} row scaling is defined as 



column scaling as C(A) = {yJ^ — j. assuming that all denominators are nonzero. The iterative 

process ...CRCR(A) is called Sinkhorn's iterative scaling (SI). In terms of the matrix Ax the 
map HS(X) can be realized as follows: 



So, the map HS(X.) is indeed a (rather far-reaching) generalization of Sinkhorn's scaling. Other 
generalizations (not all hyperbolic) can be found in [23], 0, 0. I 

Lemma 2.10 from |19| allows to use (HSI) to solve Problem 1 for P-hyperbolic polynomials q 
in the same way as it was done for the perfect matching problem in j^S] , j5U] ; and for the 
Edmonds' problem in [2] . The corresponding complexity is 0(n log(g(e))) iterations of (HSI) 
, which can be done in 0(n 3 log(g(e))) oracle calls . The algorithm works in the following way : 
Run K = 0(n \og(q(e))) Hyperbolic Sinkhorn Iterations Xj + i = HS(X-j) ; if DS(X.{) < ^ for 
some i < K then the p-mixed form M p (Xq) > , and M p (Xq) = otherwise . 



Ahsqq = C(R(Ax)) 
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5 Half-Plane Property 



The following definition is from [7j. 

Definition 5.1: A polynomial P{z\, z n ) in n complex variables is said to have the "half- 
plane property" if P(z%, z n ) ^ provided Re(zi) > . I 

In a control theory literature (see jSH] ) the same property is called Wide sense stability . And 

Strict sense stability means that 

P( Zi, Zn) provided Re(zi) > . 

The following simple fact shows that for homogeneous polynomials the "half-plane property" 
is , up to a single factor , the same as P-hyperbolicity . 

Fact 5.2: A homogeneous polynomial R(z%, z n ) has the "half-plane" property if and only if 
the exists real a such that the polynomial e ia R{z\, z n ) is P-hyperbolic polynomial with real 
nonnegative coefficients . I 

Proof: 

1. Suppose that R(zi, z n ) = e~ ia Q(zi, z n ) where a is real and Q is P-hyperbolic. 
Then Q is (1, 1, )-hyperbolic and all real vectors (x\, x n ) with positive coordi- 
nates are (1, 1, )-positive . Therefore Q is {x\, x n )-hyperbolic for all real vectors 
(xi, ...,x n ) £ R++ with positive coordinates .It follows that \R(x± + iyi, ...,x n + iy n )\ = 
\Q(xi + iyi, ...,x n + iy n )\ = \Q(xi, x n ) UiKkKni 1 + iX k)\ , where (Ai,...,A n ) are real 
roots of the real vector (yi, y n ) in the direction (xi, ...,x n ). 

This gives the following inequality , which is equivalent to the " half-plane property" of R 

\R(x! + iyi, ...,x n + iy n )\ > \R{x\, ...,x n )\ = (23) 
= \Q(xx,...,x n )\ > : 
(xi,...,x n ) e R\ + ,(y 1 ,...,y n ) € R n 

2. Suppose that R(z\, z n ) has the "half-plane property" and consider the roots of the 
following polynomial equation in one complex variable : P(x\ — z, xi — z, ...,x n — z) = , 
where (xi, x n ) £ R n is a real vector , z = x + iy S C. If the imaginery part Im{z) = y 
is not zero then , using the homogeniuty , R(i x ~^ 1 + 1, ...,i x ~* n + 1) = , which is 
impossible as R has the "half-plane property". Therefore all roots of R(X — te) = 
are real for all real vectors X £ R n (here e = (1,1, ...,1)). In the same way all roots 
of R{X — te) = are real positive numbers if X € R++ ■ It follows that if X G R n 
then R(X) = R{e) Y[l < k < nA^(A) , where (Ai, A n ) are ( real ) roots of the equation 
R(X — te) = . Thus the polynomial ("^y)-^ takes real values on R n and therefore 
its coefficients are real . In other words , the polynomial (ftj^j)R is P-hyperbolic . If 
R(l, 1, 1) = e~ ia \R(l, 1, 1)| then the polynomial e ia R is also P-hyperbolic . 
(Recall that the coefficients of any P-hyperbolic polynomial p are nonnegative for they 
are p-mixed forms of e-nonnegative tuples , and p-mixed forms of e-nonnegative tuples 
are nonnegative if p(e) > |2"fij.) 
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We use this observation to show that Theorem 2.2 in this paper implies (and seriously 
strengthens) Theorem 7.2 in [7] , which is the main result of a very long recent paper [7] . 

5.1 Submodularity and hyperbolicity 

Let p be a P-hyperbolic polynomial of degree n in n variables . It follows from Theorem 2.2 
that r = (ri, r%, r n ) € supp(p) if and only if the following inequalities hold : 

r(S) = J2 r i^ R ( S ) = Rank p (J2ei);S C {1,2, ...,n}. 

Fact 5.3: The functional R(S) = Rank p (J2 ieS ei) is normalized , i.e. i?(0) = , and submod- 
ular , i.e. R(A UB)< R(A) + R(B) - R(A D B) : A, B C {1, 2, n} . I 

Proof: Associate with two subsets A, B C {1, 2, n} the following three e-nonnegative vectors 

x = Y,ieA\(AT)B) e i> U = YsitADB e i> z = J2ieB\(AnB) e i- 

We need to prove the inequality Rank p (x + y + z) < Rank p (x) + Rank p (z) — Rank p {y). This 
inequality is obvious and well known for positive semidefinite matrices . The extension to e- 
nonnegative vectors respect to e-hyperbolic polynomial p is done in the same way as in the 
proof of Corollary A. 3 : consider a hyperbolic in the direction (1, 1, 1) polynomial 

L(cei, a 2 , a 3 ) = M p (k, k, e, e), k = a\x + a 2 y + a 3 z; 

where the vectors x,y,z are e-nonnegative respect to hyperbolic polynomial p , and the tuple 
(k, k, e, e) consists of Rank p (x + y + z) copies of k and n — Rank p (x + y + z) copies of e. 
After that apply Theorem 1.5 . I 

Corollary 5.4: 

1. A support supp{p) of P -hyperbolic polynomial p is an intersection of the integral polyma- 
troid {(ri,...,r n ) : r(S) = Ei 6 5 r * - R ( S ) = Rank p (J2ieS e i)'^ S c {M, •••,™}} with the 
hyperplane {(n,...,r n ) : J2i<i<n r i = n }- 

2. A support supp(R) of any polynomial R with the "half-plane" property is a jump system . 

Proof: (Consult |24j for a definition and some properties of jump systems and integral poly- 
matroids ) . This Corollary follows directly Theorem 2.2 , Fact 5.3 and Proposition (3.1) in 

m ■ 1 

It is quite amazing how the two communities , "hyperbolic" and "half-plane" , were not aware 
about each other results for a long , long time . (Interestingly , two authors of |7J and one 
author of j^TJ were with the same department until very recently . Perhaps , one needs to be a 
dilettante to notice a bridge .) 



17 



6 Conclusion and Acknowledgments 



Univariate polynomials with real roots appear quite often in modern combinatorics , especially 
in the context of integer polytopes . We discovered in this paper rather unexpected and very 
likely far-reaching connections between hyperbolic polynomials and many classical combinato- 
rial and algorithmic problems . (The author taught about "On hyperbolic nature of perfect 
marriages" as a title of this paper , but with the current climate it could be understood in 
many ways .) There are still several open problems . The most interesting is Conjecture 2.11 
in this paper , which is a generalization of the van der Waerden conjecture for permanents of 
doubly stochastic matrices and many others related questions . 

For a hyperbolic in direction (1, 1, .., 1) polynomial Mul(yi, y Tl ) = y\yi---y n Conjecture 2.11 
is equivalent to the famous van der Waerden conjecture for permanents of doubly stochastic 
matrices , proved in |15j , ^S] . For a hyperbolic in direction / polynomial det(X) , X is n x n 
hermitian matrix , it is equivalent to Bapat's conjecture [£] (it was also hinted in ^H] ) , proved 
by the author in , [33] . It also holds for the Moore determinant Det{M){Y) , Y is n x n 
quaternionic hermitian matrix , with the proof essentially the same as in |35| . 
Another , equivalent form of "hyperbolic" (or "half-plane" ) van der Waerden conjecture can 
be formulated as follows : 

Conjecture 6.1: Consider a homogeneous polynomial p(z\, z n ) of degree n in n complex 
variables . Assume that this polynomial satisfies the property : 

\p(zi,...,z n )\ > Hi<i< n Re(zi) on the domain {(zi,...,z n ) : Re(zi) > 0, 1 < i < n} . 

Is it true that I „ 8 " p \ > i*L ? 
(Notice that Theorem 2.10 and Fact 5.2 imply that dz ^ dz J > / .) I 

It would be very interesting and enlighting to prove Conjecture 2.11 using methods of the 
theory of functions of many complex variables. Fact 5.2 , together with other results in this 
paper , makes a connection between the Complexity Theory and the theory of linear time- 
miltidimensional systems : all "hard" instances of Problem 1 are necessary unstable polynomi- 
als. 

Another interesting conjecture is related to the majorization : 

Conjecture 6.2: Consider the doubly-stochastic and P- hyperbolic homogeneous polynomial 
p(x\, ...,x n ) of degree n in n real variables . 

Let A(A) G R n be a real n-dimensional vector , whose coordinates are the roots of the equation 
p(X — te) = 0, where X E R n and e is the vector of all ones . Then there exists a n x n doubly 
stochastic matrix A such that A(X) = AX. 

(Some partial and related results in this direction can be found in ; this conjecture is 
true for determinantal polynomials .) I 
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A natural extension of Problem 1 is for P- hyperbolic polynomials to approximate g Xl g x ^ -p(xi,. 
within a multiplicative factor using deterministic ( or randomized ) polynomial-time oracle al- 
gorithms . It is not clear to the author whether known recent randomized algorithms for (1 + e) 
approximation of the permanent Per(B) of entry-wise nonnegative matrix B can be done in 
the "oracle fashion" , i.e. using only some outputs of the multilinear polynomial 

q(x 1 ,...,x n ) = Yl B (^j) x j- 

l<i<n l<j<n 

If hyperbolic van der Waerden conjecture is true then the technique in this paper , similarly to 
|20| and |2j , [22] , would produce a deterministic polynomial-time oracle algorithm with ^ 
multiplicative factor . 

The technique developed in this paper can be applied to other "noble" desicion problems 
. For instance , checking factorizability of P-hyperbolic polynomials can be also done in de- 
terministic oracle polynomial time . The factorizability is closely related to the hyperbolic 
generalization of the indecomposability of matrix tuples |22| . 

This paper is probably the first one which uses Theorem 1.5 in the combinatorial context . We 
expect many more such applications of Theorem 1.5 . This (very nontrivial) theorem , when in 
good hands , is a powerful tool allowing reasonably simple and short proofs . 

I would like to acknowledge a great influence of amazingly clear paper |2H] . It is my pleasure 
to thank Adrian Lewis for numerous as e-mail as well phone communications. Many thanks to 
the fantastic library of Los Alamos National Laboratory and Google. 
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A Proof of the (main ) Theorem 2.2 

Before proving Theorem 2.2 , we will recall some basic properties of p-mixed forms and prove 
a few auxiliary results . The following fact was proved in [2Sj 
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Fact A.l: Consider a homogeneous polynomial p(x),x £ R m of degree n in m real variables 
which is hyperbolic in the direction e. Then the following properties hold . 

1. The p-mixed form M p (x\, .., x n ) is linear in each Xj, 1 < i < n. 

2. If xi, X2, .., x n -i are e-nonnegative then the linear functional l{x) = M p (x±, .., x n -±, x) is 
nonnegative on the closed cone N e of e-nonnegative vectors . 

3. If the tuples (x±, ..,x n ), (yi, -.,2/™), (xi — y±, ..,x n — y n ) are e-nonnegative then 

< M p (y u ..,y n ) < M p (x!, ..,x n ). 

4. Fix e-positive vector d and consider the following homogeneous polynomial Pd(x), x £ i? m 
of degree n — 1 in m real variables : Pd(x) =: M p (x, x, x, d) . Then is hyperbolic in 
any e-positive direction v 6 C e (p) . If g 6 C^p) ( e-positive respect to the polynomial p 
) then also g G C v (pd) for all u G C e (p) . 

I 

The next fact is well known . 

Fact A. 2: Consider a sequence of univariate polynomials of the same degree n : Pk(t) = 
J2o<i<n a iM l ■ suppose that linifc^oo a i>k = ^,0 < i < n and a n ^ . 

Define P(t) = J2o<i<n a i^ 1 ■ Then roots of converge to roots of P . In particular if roots of 
all polynomials Pk are real then also roots of P are real ; if roots of all polynomials P^ are real 
nonnegative numbers then also also roots of P are real nonnegative numbers . I 

The following corollary of Theorem 1.5 plays crucial role in our proof of Theorem 2.2 . 

Corollary A. 3: 

1. Consider a homogeneous polynomial p(x),x £ R m of degree n in m real variables which 
is hyperbolic in the direction e. Let x\,X2,x% be three e-nonnegative vectors and d = 
x\ + X2 + X3 is e-positive . Assume wlog that p{x\ + X2 + £3) = 1 • Then there exists 
three symmetric positive semidefinite matrices A, B, C such that p(a±xi + 02X2 + 03^3) = 
det(aiA + CI2B + CI3CO for all real 0,1,0,2,0,3. Additionally , the roots 0/01X1 + 02X2 + 03X3 
in the direction d , i.e. the roots of the equation p(a\X\ + 02X2 + 03X3 — td) = , coincide 
with the eigenvalues of a±A + 02-B + a%C . 

2. Theorem 2.2 is true for e-nonnegative tuples 

(X) = (xi, ...x n ),Xi £ R m consisting of at most three distinct components , i.e the cardi- 
nality of the set {x\, ...x n } is at most three . 

Proof: 
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1. Consider the following homogeneous polynomial L(bi,b 2 ,b 3 ) = p(bix\ + 62^2 + b 3 (xi + 
X2+x 3 )) °f degree n in 3 real variables . It follows from Theorem 1.5 that there exists two 
real symmetric matrices A and B such that L(b±, 62, ^3) = det(&i A+&2-B + &3-O . It follows 
that they both positive semidefinite , and C = I — A — B is also positive semidefinite . 
Take a real linear combination z = a\X\ + 02X2 + 03X3. Then 

p(z - t(xi + X 2 + X 3 )) = 

det((ai - a 3 )A + (a 2 - a 3 )B + a 3 I - tl) = 
= det(oi^ + a 2 B + a 3 C - tl). 

This proves that p(a±xi + 0,2X2 + 03^3) = det(a\A + 02-6 + 03c) for all real a±, 0,2, a 3 by 
putting t = 0. And it also proves the "eigenvalues " statement . 

2. Consider e-nonnegative tuple (X) consisting of rj copies of X{ , 1 < i < 3 ; r\ + T2 + r 3 = n 
. Assume that d = x\ + X2 + x 3 is e-positive (if it is not then M p (X) = by a simple 
argument based on the monotonicity of p- mixed forms ) . It follows from the polarization 
formula (10) that 

M p (X) = dip(t 1 jX 1 +t 2 ,iX2+t3 ! iXs), 
l<j<fc<oo 

and this formula is universal , i.e. holds for all homogeneous polynomial of degree n , in 
particular for det(A) , X is n x n symmetric matrix . Therefore , using the first part 
of this Corollary we get that the p-mixed form M p (X) = D(A) , where the matrix tuple 
A consists of r\ copies of A , T2 copies of B and r 3 copies of C and D(A) is the mixed 
discriminant . Using Rado theorem for mixed discriminants we get that D(A) > iff 

RankQ^Ai) > for all Sc {1,2,3}. 

i&S i&S 

But from the first part we get that Rank(^2 i&s Ai) is equal to p-rank Rank p (J2i<=s x i) 
ofV, s .r, lor all ,V : {1.2.3} . 

I 

Proposition A. 4: Consider similarly to part 4 of Fact A.l the polynomial p^ix) =: M p (x, x, x, 
where d is e-nonnegative and Rank p (d) > 1 . Then pd is hyperbolic in any direction z G N e (p) 
which is e-nonnegative and satisfies the following inequalities : 

Rankp(z) > n — 1; Rank p (z + d) = n. (24) 

Also , if y £ N e (p) is e-nonnegative then also y £ N z {jpd) , i.e. is z-nonnegative respect to the 
polynomial pd- 

Proof: Let z £ N e {p) be e-nonnegative vector satisfying (24) . Consider univariate polynomial 
P(t) = M p (tz+x, tz+x, ...,tz+x, d) . Then P(t) = J2o<i<n-i an d o-n-i = M p (z, z, z, d) . 
It follows from Corollary A. 3 that a n -\ > . Consider now a sequence of univariate polynomials 
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Pk(t) = M p (tzk + x,tzk + x,...,tzk + x,dk) ■ Where Zk,dk are e-positive and lim^oo zt = 
z, linifc^oo dk = d . Then the coefficients of polynomials converge to the coefficients of the 
polynomial P . It follows from part 4 of Fact A.l that the roots of Pk are real . Since a n _i > 
hence using Fact A. 2 we get that the roots of P are also real . This exactly means that the 
polynomial pd is hyperbolic in direction z . The d-nonnegativity statement follows from the 
nonnegativity part of Fact A. 2 . | 

We are ready now to present our proof of Theorem 2.2 . The proof is by induction in the 
degree n . The main trick which we used is that to justify the induction , i.e. that if the 
generalized Rado conditions hold for hyperbolic polynomial p of degree n then the generalized 
Rado conditions hold for some auxiliary hyperbolic polynomial pd of degree n — 1 , we need 
to prove Theorem 2.2 for tuples consisting of at most three distinct components . And this 
particular case follows from the classical Rado theorem via Theorem 1.5 and Corollary A. 3 . 

Proof: (Proof of Theorem 2.2 ) . 

The "only if part is simple . Indeed supposed that there exists a subset S C {l,2,...,n} 
such that Rank p (J2i£s x i) < 1^1 > i- e - usm g the identities (14) M p (k,k, ...k,d, ..,d) = , where 
k = J2ies x i > d ^ Ce(p) is e-positive and the n-tuple (k,k, ...k,d, ..,d) consists of \S\ copies 
of k = J2ies x i ■ Let d be any e-positive positive vector such that d — X{ is e-nonnegative , 
1 < % < n . Using the monotonicity of p-mixed forms we get that 

M p (x\, x n ) < M p (k, k, ...k, d, .., d) = 0. 

Our proof of the "if part is by induction in the degree n . Suppose that the generalized Rado 
conditions (15) hold . Then at least Rank p (x n ) > 1 . Consider the following homogeneous 
polynomial of degree n — 1 : 

p A (x) = M p (x,x, ...,x,d), d = x n . 
We get from Proposition A. 4 the following assertion : 

The polynomial pd{x) is hyperbolic in direction z = Ei<Kn-i x i an d the vectors X{ G N z (pd), 1 < 
i < n — 1 , i.e. are z-nonnegative respect to the polynomial pd- 

Indeed , it follows from the generalized Rado conditions (15) that Rank p (z) > n — 1 and 
Rank p (z + d) = Rank p (J2i<i< n xi) = n . 

Next we show that the n— 1-tuple Y = (x±, ...,x n -i) satisfies the generalized Rado conditions 
for z-hyperbolic polynomial pd of degree n — 1 : 

Rank Pd (J2 x i)>\S\ for all S C {1, 2, n - 1}. 

i&S 

Or equivalently (see formulas (14) ) , that 

M p (k, .., k, z, z, d) > 0; k = J^ies x ^ ( 25 ) 
z = Ei<k„-i Xi, d = x n ,S C {1, n - 1}, (26) 

where the n-tuple T = (k, .., k, z, z, d) consists of \S\ copies of k , n — 1 — \S\ copies of z and 
one copy of d . 
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It is easy to see that the generalized Rado conditions for the n-tuple T are implied by the 
generalized Rado conditions for the original n-tuple X = (xi, ...,x n -i,x n ) . Since the n-tuple 
(k, .., k, z, z, d) consists of at most three distinct components hence we can apply part 2 of 
Corollary A. 3 . Therefore we get that indeed 

M p (k, .., k, z, z,d) > and hence the following inequalities hold : 

Rank Pd (^Xi) > \S\ for all S C {1, 2, n - 1}. (27) 

ieS 

Thus , by induction in the degree , we get that pj- mixed form M Pd (xi, ...,x n -±) > : 
the polynomial pd of degree n — 1 in m real variables is z-hyperbolic . But 

M Pd (xi,...,x n _i) = dai dn d ^ n l Pd(Ei<i<n^i atiXi) = 

= <9cti...c9a n _i ^?El<Kn-l Oi{Xi, El<i<n-1 O-i^i, X n ) = (n — l)\M p (xi, ...,X n ). 

We conclude that if Theorem 2.2 is true for n — 1 then it is also true for n , and the case 
"n = 1" is trivially true . | 

Remark A. 5: Consider a mixed discriminant D(A) , where A = (Ai,...,A n ) is a n-tuple of 
positive semidefinite n x n hermitian matrices , i.e. A{ y . Recall that in this case -D(A) > 
; and D(A) > if and only if there exists n linearly independent vectors vi,...,v n such that 
Vi 6 Im(Ai), 1 < i < n. 

In the proof of Theorem 2.2 we encountered the following tuple of positive semidefinite matrices 

A = (A, A, B, B, C) consisting of I copies of A , m copies of B and one copy of C . 
Moreover , this tuple is even more special . I.e. B — A y , Rank(B) > n — 1 , Rank(A) > I , 
rank(C) > 1 ,Rank(A + C) > I + 1 and Rank(B + C) = n . 

For such tuples the Rado theorem has very elementary proof , which we sketch below . 

There are two cases . First case is when Rank(B) = n , it is simple and left to the reader . 
Second case is when Rank(B) = n — 1 . 

This is how in this case we can choose vectors v n 6 Im(C);vi, ...,vi G Im(A);vi + i, ...,v n -\ G 
Im(B) in such a way that (v±, v n ) is a basis : first choose nonzero v n £ Im{C) which does 
not belong to Im(B) , second choose any I linearly independent vectors v±,...,vi G Im(A , 
third choose any n — l — 1 linearly independent vectors in Im{B) n L(v±, vi) 1 - . (L(v±, vi) 1 - 
is a linear subspacespace orthogonal to the linear subspace L(v±, vi) which is spanned by 
(vi, -,vi).) 
I 

B Proof of Proposition 2.6 

Proof: Assume wlog that q{a\, ...,a n ) = 1 . It follows from the Euler's identity that 

^2 ai-^— <?(«!,•••, «n) = n. 
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Let q(a>i, ...,a n ) — J2( ri ,...,r„)esu PP (q) a (n,...,r n ) Ui<i< n a 7 ■ 
Define the following nonnegative real numbers : 

&(n,..,r„) = a (n,...,r n ) ]J a? ,(n, ...,r n ) € supp(q). 

l<i<n 

Then a i ^-q( y a 1 ,...,a n ) = E(r u ...,r n )esupp( q )rib( ri ,...,r n ) ■ 

Suppose that for some subset S C {1, 2, n}, 1 < \S\ < n we have the inequality J2ies r, i < 
jS*! for all (ri,...,r n ) £ supp(q) . Then^2 ieS ai-£^q(cti, ...,a n ) < | *S'| — 1 . But the condition (17) 
says that ai-£^-q(ai, ...,a n ) = 1 + 5$ and Xa<i< ra l^i| 2 ^ ^ • By the Cauchy-Schwarz inequality 

, E ie5 N < Vi? < 1 • Therefore, 

XI ai—q(ai, a„) > |S| - ^ | <5» | > \S\ - 1. 
The last inequality gives a contradiction . | 



C A sketch of a proof of Corollary 2.5 

Proof: By Theorem 2.2 the conditions (1) and (2) are equivalent . (2) implies (3) for any 
homogeneous polynomial with nonnegative coefficients . 

Let oti = e Vl , 1 < i < n; J2i<i<n Hi = 0- Consider the following convex functional 

f( yi ,...,y n ) = log(q(e^,ey\...,e^)). 

Here q(x),x G R n is a homogeneous polynomial of degree n in n real variables with nonnegative 
coefficients . Then 



ai^q{ai,...,a n ) Q 



f(yi, ■■;Vn), l<i<n. 



q(a!,...,a n ) % 
Notice the condition (3) is equivalent to the following condition : 

hif /(yi,...,y n ) =L> -oo. 

J/l+...+J/n=0 

Consider the anti-gradient flow , i.e. the system of differential equations 
Vi{t)' = -.fn) " 1), J/<(0) = 0; 1 < i < n. 

It is well known that in this convex case the gradient flow is defined for all t > . Using the 
Euler's identity we get that 

1 -f{yi{t),...,y n (t)) = -p(t)=:- 2^ | ^ 1 

dt ^ q(a!,...,a n ) 
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It is easy to see that , because of the convexity of / , a nonnegative function (5{t) is non- 
increasing on [0, oo) . 

As infj, 1+ ... + j /n=0 /(yi,...,y n ) = L > -oo thus J °° (3{t)dt < oo . Thus lim^ 00 /3(t) = . This 
proves the implication (3) — > (4) for all homogeneous polynomials of degree n in n real variables 
with nonnegative coefficients . 

The implication (4) — ► (5) is obvious . The implication (5) — ► (6) for general homogeneous 
polynomials of degree n in n real variables with nonnegative coefficients is Proposition 2.6 . 
Finally , the implication (6) — ► (2) follows fairly directly from Theorem 2.2 . | 



D Lower bounds on the number of oracle calls for the exact 

dxi-.dxn 1 



computation of «A-p(ii, ...,x n ) 



Definition D.l: Call a set {X\, ...,X m },Xi G C n e-universal if there exist complex numbers 
ci, c m such that for any homogeneous polynomial p(.) of degree n in n complex variables the 
following inequality holds 

I dxZdx n P( X l> X ") ~ T,l<i<rn CiP(Xi)\ (28) 

< emax( ri) ... !rn )ei n ,n W u ...,r k \, 

where a, ri ,...,rni 

(n, ...,r n ) e I n>n are the coefficients of the polynomial p(.) . | 
Lemma D.2: // the set {Xi, ...,X m },Xi G C n is O-universal then 

m -W^W"^M (29) 

If the set {Xi, X m }, Xi £ C n is e-universal then 

1 n' 

m ^ min([ 7 ] ' [f]!(n-[§])! ) (30) 

Proof: Define a monomial M ri rjl (ari, x n ) As {Xi,...,A m } is universal 

thus the exists complex numbers (ci, c m ) , which are wlog are all nonzero , such that 

Y, CtMr^M Xi) = 
l<i<m 

if (ri,...,r n ) G /(n, n), (n, r n ) / (1,1,.. .,1) ; 

i A i 

and Z)i<i<m Ml,i,...,i(c" -Xi) = 1 ; define Yi = c" Xj (here c" is one of the nth complex roots of 

Q )• 

Let Half = {(ri,...,r n ) : ti € {0, 1}, J2i<i< n r i = [§] • Notice that the cardinality \Half\ = 

[^.(n-m — K ■ 

Define the following two K x m complex matrices : 

W((r 1 ,...,r n ),j)=M ri _ rn (Y j ), 
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V({n,...,r n ),j) =M 1 _ rii ..., 1 _ rn (y i ) : 
(n, •••,r n ) e flaZ/, 1 < J < m. 

Clearly , i?an£;(W) = Rankiy) < m . On the other hand the O-universality condition implies 
the matrix identity 

WV T = I (31) 
Therefore m > Rank(W) > \Half\ = ^n]\(n-[^])\ ■ 

If the set {Xi, X m } is e-universal and de < l,d G N then Rank(WV T ) > d . This proves 
(30) . | 

Remark D.3: The identity (7) is a particular case of a slightly more general one : 

— — ^—p(xi,...,x n ) = E(p(z 1 ,z 2 ,...,z n ) Y[ Zi), ( 32 ) 

l<i<n 

where (z±, z 2 , z n ) are independent complex random variables such that E(zi) = and 
E{z{Zi) = 1 for all 1 < i < n . The identity is easily proved by checking it for all mono- 
mials M rit ,„ t r n , (n, ...,r n ) £ I(n,n) . If p(.) is a multilinear polynomial , i.e. 

p(x 1 ,...,x n ) = Yl ( E A (hj) x j), 

l<i<n l<j<n 

then Q Xl d "d x p(%i, x n ) = Per(A) , where per(A) is the permanent of the matrix A . Clearly 
, lower bound m > p j j^-if" ] )! also holds for multilinear polynomials and even for powers 

(12i<i<n a i x i) n ■ It is very likely that the actual lower bound is 2 n_1 and that it does exist 
somewhere in the geometrical designs literature . In the case of permanents , the formula (7) is 
essentially the Ryser's formula [3]; and Lemma D.2 says that , in some sense , it is an optimal 
formula for computing permanents . 

Another equivalent formulation of Lemma D.2 is the following statement : 
Let a set of complex vectors 

S = {X l = (x lA , ...,x l>n ) eC n :l<l< . (2n ~^ )! . } 

(n — l)!n! 

be a Haar set for the monomials M rir __ jrn : (n, ...,r n ) £ 

I.e. the square matrix {M ri) ... irri (Xj) : X{ 6 S; (n, ...,r n ) 6 I n ,n\ is nonsigular . 
If 

[] x hi = J2 c i( E Y(k,i) Xi ) n 

l<i<n l<k<m l<i<n 

for all 1 < I < fe^i and some complex numbers {c/%; Y(k, i) : 1 < k < m, 1 < i < n} then 

"[S]!(n-[S])! y/n' 
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